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The interplay between quarter-filled and half-filled umklapp scattering has been examined by 
applying the renormalization group method to a one-dimensional quarter-filled electron system 
with dimerization, on-site (U) and nearest-neighbor (V) repulsive interactions. The phase dia- 
gram on the U-V plane is obtained at absolute zero temperature where the Mott insulator (the 
charge ordered insulator) is found for smaller (larger) V. By choosing the moderate parameter 
in the region of Mott insulator, it is shown that the resistivity exhibits a crossover from behavior 
of quarter-filling to that of half-filling with decreasing temperature. 
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The organic conductors, Bechgaard salts, which 
show a crystal structure with a stacking of organic 
molecules along one-dimensional axis, form a quasi-on 
dimensional system with a quarter-filled hole bando 
Recently, the states above the phase transition tempera- 
ture have been studied extensively to show several exper- 
imental observations which may be relevant to the uncon- 
ventional non-Fermi liquid. Optical and photoemission 
experiments exhibit_hhe charge gap in both TMTTF salts 
and TMTSF salts-Mlf The charge ordered (CO) insu- 
lating state in (TMTTF) 2 AsF 6 and (TMTTE) 2 PF 6 has 
been maintained by the NMR measuremeni£l ) and the 
huge dielectric response in (TMTTF) 2 PF 6 B' indicates 
also the CO state. It is necessary to understand these 
experimental facts in consistent with a general phase di- 
agram of ihe Bechgaard salts indicating a dimensional 
crossover .EP 

Two mechanisms Jjpr the insulating state has been pro- 
posed theoretically!? One is of them is the-half-filled 
umklapp scatterings induced by dimerizationtf and the 
charge gap has been examined for the one-dimensior±al 
quarter-filled Hubbard model with the dimerization.Ef 
The other is the commcnsurability at the quarter-filling. 
The quarter-filled extended Hubbard model without the 
dimerization has been examined by using a method of nu- 
merical diagonalization, where the phase diagram shows 
the insulating state., with the charge gap for a large 
Coulomb repulsion. c3i3 The similar phase diagram has 
been obtained analytically by calculating the magnitude 
of the umklapp scattering at quarter-filling and by us- 
ing the bosonizatiatL method and the renormalization 
group technique.t 3 Jll-3' The model treating both of these 
two umklapp scattering on the same footing is required 
for investigating Bechgaard salts. By using the DMRG 
method, the charge gap of the model has been estimated 
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where the metallic state is reduced to the insulator in 
the presence of dimerization.E^P 

In the present letter, we investigate the interplay of the 
half-filled umklapp scattering andn-the quarter-filled one 
by extending the previous workc3llj' to the case includ- 
ing the dimerization. It is demonstrated that the com- 
petition between these two scatterings gives rise to two- 
kinds of unconventional insulating states in the ground 
state and that a crossover between these two states ap- 
pears in the resistivity at finite temperatures. 

We consider a one-dimensional quarter-filled extended 
Hubbard model with dimerization. The Hamiltonian is 
given by 



H = - 



j2(t+(-iyt d ) (ct a c i+1)(T + h.c. 



' X! "i.T n j,l + V Y1 n i n i+ l 



(1) 



where n 



cJ )CT Cj i<7 and n 



3 = n 0A + n 3,l- 



The fcrmion 



operator ct a denotes a creation of the electron at the j-th 
site with spin a. We have two kinds of transfer integrals 
t+t<i and t— <d due to the dimerization. The quantities U 
and V denote the on-site and nearest-neighbor Coulomb 
repulsive interaction. In the presence of the dimerization, 
the unit cell is given by two lattice sites with 2a where 
a is the lattice constant. By introducing operators a 
and cjj, a a which are those for the even and odd sites in 
the n-th unit cell (R n = 2an), respectively, we perform 

-ikR n „A 

, where Nq = 



the Fourier transform c^ c 



N, 



-1/2 



and C* = N 



- AT- 1 / 2 



j: No e- 



ik(R n 



J R n +a,a 



L/(2a) and L is the length of the system. Then the 
kinetic term of the Hamiltonian, TCq, is diagonalized as 
= J2k.cr £ k [d{ M d k ,a - ul a u k ,<j] with a dispersion re- 
lation, Ek = — 2[t 2 cos 2 ka + t 2 sin 2 ka] 1 / 2 . In a diago- 
nalized basis, the lower and upper band fermion opera- 



tors are given by dk.o = [e 



]/V2 and 
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Fig. 1. Schematic view of the two limiting cases for the insulating 
state: the Mott insulating state (a) and the CO insulating state 
(b). The dashed wave indicates the 4fep charge density wave 
with Op = tt/4 (a) and P = or tt/2 (b). 



Uk,a = [e l9fe a — e~ iek cf a ]/V2, respectively, where 
tan20fe = — (td/t) t&nka. Since a gap 4td appears at 
k = ±ir/2a and the first Brillouin zone is reduced to half, 
the system can be regarded effectively as half-filling with 
fcp = 7r/4a. 

In order to obtain the effective Hamiltonian for the 
states near ±fcp, one integrates out the contribution 
from the upper band, which leads to the quarter- 
filled umklapp scattering.li3 We apply the bosoniza- 
tion method to electrons in the lower band with the dis- 
persion linearized as £k — > ±Vp(\k\ — kp) where vp = 
V2ta[l - (t d /t) 2 ]/[l + {td/t) 2 } 1 / 2 . Then the Hamilto- 
nian, which is rewritten in terms of the charge and spin 
phase variable, is decoupled into two parts of the charge 
and spin degr. 
expressed asli^ 



;es of freedom. The charge part, Ti. p 



is 



(d x 9 p f + K p (d x <f> p ) 2 

K p 



2^ / <Lr " in2 ^' 



+ 



9l/4 



dx cos 49 n 



(2) 



where a is a cutoff of the order of lattice constant and 
[0 p (x), (j} p {x')} = i7rsgn(a; — a;'). In eq. (2), the charge 
velocity and the Tomonaga-Luttinger exponent are given 
by v p = v F [(l + g 4p /2TTv F ) 2 - (g p /2Trvp) 2 ] 1/2 and K p = 
[(2irv F +g ip -g p )/ (2Trvp+g ip +g p )} 1 / 2 . The magnitude of 
half-filled umklapp scattering, 51/2, and that of quarter- 
filled one, 31/4, are given by 

Ua-^—U(U-2V)I 2 (A) 
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where A ee [1- (t d /i) 2 ]/[l + (t d /f) 2 ] and B = [2t d /t]/[l- 
(tfjlt) 2 ]. The quantity g^u has been derived for td 
0.113 The constants g p and g\ p are given byE 



V 



9p 



04p 



(U + W)a 
Aa 2 - 



4irv F 
(U + AV)a 



(u - 2vy + u z + (2vy i x {a) 



Aa 2 

4:TtV F 



(U - 2V) 2 + U 2 B 2 + (2V) 2 ] I 2 (A) 



where h(A) = Io /2 ^Y^=±[2 
and I 2 {A) ee 2/ V2 d(^ [2- 



2^1" 



eA cos ip\ 1 , 
Acos v?] _1 - 



■ „ - , n . r _ rVl + A COS (f+y/l 

These quantities are calculated as h(A) ~ ii(l) = 
ln(V2 + 1) and I 2 (A) ~ I 2 (l) = 2(^/2 + 1). The or- 
der parameters p£ the 4fcp charge density wave (CDW) 
state is given h$3 



O 



4fc F -CDW 



cx cos (Ak F x + 29 p ) 



(5) 



From eq. (2), one finds the following facts. In the case 



of 



0, 9 D is locked to o 



7r/4 for gi/2 > 0, 



(or 9 P = 37r/4 for g x j 2 < 0) with < 9 p < n. In 
this case, the charge configuration for 4fcp-component 
is shown schematically in Fig. 1(a), which corresponds 
to the Mott insulating state. We note that the sign of 
the half-filled umklapp scattering depends on that of the 
dimerization, i.e., gi/ 2 > for td > and g±/ 2 < for 
td < 0. On the other hand, in the case of gi/ 2 = 0, the 
phase is locked to 9 p — 0, n/2 for g 1 i i < or 9 p = tt/4, 
37r/4 for g 1 / i > 0. The locking of 9 p for g 1 / i < (i.e., 
U < 4V) corresponds to the CO insulating state shown 
in Fig. 1(b). In the presence of both gx/ 2 and <?i/4, the 
behavior becomes quite different depending on the sign 
of <?i/4. In the case, g 1 / /l > 0, the two kinds of umklapp 
scattering are compatible and g^-term removes the de- 
generacy between 9 p = tt/4 and 9 P = 37r/4. However, in 
3i/4 < case, the two kinds of umklapp scattering com- 
pete with each other since the locking positions of 9 p 
are different for respective umklapp scattering. The CO 
state, which appears in large V, has been obtained by 
the mean-field approximation methodlilP and the com- 
petition between these two kinds oLpotentials has been 
shown in terms of phase variableJlS* Here we investi- 
gate the competition between the Mott insulator and the 
CO insulator, by employing the renormalization group 
method, where the quantum fluctuation is fully taken 
into account. 

The renormalization group equations for the coupling 
constants in the Hamiltonian, eq. (2), are given by 

d K n = -2G 2 



dl 



K 2 

1/2 A P 



' ^1/4 Kp 



-jy Gx/2 = (2 - 2K P ) Gx/2 
— Gx/i = (2 - 8^) Gx/4 



Gx/2 Gx/i 

X G 2 
2 G i/2 . 



(6) 
(7) 
(8) 



where G z = g z /(2irv p ) (z = 1/2, 1/4) and I = ln(a'/a) 
with the new length scale a' being larger than a. We 
note that these renormalization group equations can be 
derived in a way similar to Ref. [19]. Equations (7) and (8) 
indicate the fact that Gi/ 2 and G x /x are suppressed by 
each other for G1/4 < and are enhanced for Gx/a > 
as discussed in the last paragraph. 

The competition for gx/x < (i.e., U < 4V) is examined 
numerically. The parameter of the dimerization is taken 
as td/t =0.001 to show results clearly, although such a 
choice is rather small compared with that of the organic 
conductors (td/t ~ 0.1). In Fig. 1, the ^-dependences of 
K p (l), Gx/ 2 (l) and G x /x(l) are shown for U/t = 6 and 
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Fig. 2. The scaling flows of the coupling constants for V/t = 4.60 
(a) and 4.75 (b) where U/t = 6 and t d /t = 0.001. These are the 
typical scaling flows for the Mott insulator (a) and for the CO 
insulator (b), respectively. 



t d /t = 0.001 with fixed V/t = 4.60 (a) and 4.75 (b). For 
both cases, the exponent K p decreases monotonically in- 
dicating the insulating state. In Fig. 1(a), the half-filled 
umklapp scattering G\/ 2 {1) increases and becomes rele- 
vant leading to the Mott insulator. The sign of G\/±{1) 
changes from negative to positive due to the relevant 
Gi/ 2 (l). Figure 1(b) exhibits G\/ 2 {1) and Gi/ 4 (Z), which 
are quite different compared with those of Fig. 1(a). In 
this case, G\u{l) is relevant while G\/ 2 (l) decreases af- 
ter taking a maximum and becomes irrelevant. The rel- 
evance of G1/4 with negative value denotes the locking 
of 9 p being or ir/2 and leading to the CO insulating 
state. 

Based on Fig. 2, the phase diagram of the ground state 
is shown in Fig. 3 with fixed td/t = 0.001. In the Mott 
insulator, the phase 9 P is locked at 9 P — 7r/4, while in the 
CO insulator the phase is fixed at 9 P = or tt/2. The 
dashed curve denotes the boundary in the limit of id — * 0, 
which is equal to that of the metal-insulator transition 
in the quarter-filled extended Hubbard model.E3 In the 
presence of id, the metallic region for id = changes into 
the Mott insulating state and the CO insulating region 
is suppressed. At the boundary, the charge configura- 
tion changes from uniform (Fig. 1(a)) to CO (Fig. 1(b)) 
with increasing V. We note that the state at the bound- 
ary becomes metallic since this case corresponds to an 
Ising transition found in the sine-Gordon system with 
two kinds of noii=liii£ar terms, i.e., the gap collapses on 
the criticalityMB-B) 

Now we examine the competition at finite tempera- 
tures. We calculate the temperature dependence aLt 
resistivity using the memory functional approaches t 
which has been applied to the one-dimensional sys- 
tems by Giamarchi.Ej' The conductivity is given by 
(t(u>) = i[x(w) — x(0)]/w, where \ is the retarded 
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Fig. 3. The phase diagram of the Mott insulator and the CO 
insulator on the plane of U/t and V/t with fixed t d /t = 0.001. 
The dashed curve indicates the boundary in the limit of t d — » 0. 



current-current correlation function: x(w) = ({j',j)) u = 
-i J dx J™ dt{\j(x,t), j(0,0)]) e^*"" 5 *. The current op- 
erator j(x,t) is given by j(x,t) = —dt9+{x,t)/ir = 
—1VpK p H(x, t) where IT = —d x 9_/ (2ir). From the mem- 
ory function approach, u(lu) is rewritten asES* 



cr(w) 



2v p K p 



M{lo) 



(9) 



where M(u>) is the memory function defined by M (w) = 
wxH/[x(0) — x( w )]- The memory function is cal- 
culated perturbatively as M(u>) ~ — [((F; F)) — 
((F;F))l =0 }/[u;x(0)} where ((F;F))° U stands for the re- 
tarded correlation function in the absence of umklapp 
scattering and F = [j,Tt]. In terms of M(oj), the resis- 
tivity is expressed as 



p(T) = lim a-^Lu) = 



i2v p K p 



lim M(w) 



Equation (10) indicates the power-law behavior, 



p(T)(X J2 G l/2nT 
n=l,2 



(10) 



(11) 



In order to obtain the precise behavior of the resistivity, 
we use the solutions of the renormalization group equa- 
tions, in which the umklapp scattering becomes relevant 
at low temperature. By putting e~ l = 2%T/W with 
W = v p a~ 



P(T) 



the resistivity is obtained as 
4n 2 7r 



E 

n=l,2 



G? /2n (0 cos 2 {n 2 *K p (l)) 

x B 2 (n 2 K p (l), 1 - 2n 2 K p (l)) e~ l . (12) 

where B{x,y) is the beta function given by B(x 1 y) = 
T(x)T(y)/T(x + y). In Fig. 4, the temperature de- 
pendence of the resistivity is shown for U/t — 6 and 
td/t = 0.001 with fixed V/t = 4.60 (a) and 4.75 (b). 
The curve (a) shows the behavior of the Mott insula- 
tor, while the curve (b) represents that of the CO in- 
sulator. For the CO insulator (curve (b)), the resistiv- 
ity is determined by the quarter-filled umklapp scatter- 
ing, where the resistivity exhibits the power-law depen- 
dence given by p(T) cx T 16Kp ~ 3 around the temperature 
T/W^>1 x HT 1 . The temperature T/W ~ 5 x 10~ 2 , 
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Fig. 4. The temperature dependence of the resistivity for U/t = 6 
and t d /t = O.OOlwith fixed V/t = 4.60 (a) and 4.75 (b) where 
W denotes the cutoff of the order of the bandwidth. The curve 
(a) (curve (b)) corresponds to the behavior of the Mott insulator 
(the CO insulator). 



corresponding to a minimum of the resistivity, is of the 
order of the charge gap A p below which the resistivity- 
increases rapidly. It is expected—that the resistivity for 
T -c A p behaves as exp(A p /T)J$ although the pertur- 
bative renormalization group approach breaks down for 
T <C A p . For the case of the Mott insulator (curve 
(a)), we found a noticeable behavior due to the rele- 
vance of the half-filled umklapp scattering. There ap- 
pears a mid temperature region around T/W ~ 1 x 1CP 2 , 
where the resistivity shows power-law behavior of half- 
filling given by piT) oc T 4Kp ~ 3 . At higher temperatures 
(T/W ^,4 x 10~ 2 ), the contribution of the quarter-filled 
umklapp scattering becomes dominant at which the re- 
sistivity is similar to the case of curve (b). Such a 
crossover in the region of CO insulator can be found by 
the moderate choice of parameter with |G r 1 / 4 1 > Gi/ 2 , 
i.e., close to the boundary in Fig. 3. The similar result 
is found in the frequency dependence of the conductiv- 
ity at T = 0. While the optical conductivity of the CO 
insulating phase corresponding to curve (b) shows the 
power-law dependence with u(oj) oc cj 16A "p~ 5 above the 
frequency of the charge gap, the conductivity of the Mott 
insulator corresponding to curve (a) shows a crossover 
from cr(oj) oc ui 16Kp ~ 5 to <j(u)) oc w 4 ^ -5 as decreas- 
ing frequency. Thus, the crossover from quarter-filling 
to half-filling is found with decreasing temperature or 
frequency when parameters are located in the Mott in- 
sulating region and close to the boundary. 

Finally we comment on observation in the Bechgaard 
salts, (TMTTF) 2 X and (TMTSF) 2 X. The optical mea- 
surement in the normal states of these salts shows that 
the optical conductivity has a power-law dependence as 
<j(oS) oc a;™ 1,3 at high frequency, indicating a fact that 
the quarter-filled umklapp scattering is dominant .Ef 
Our calculation suggests that the quarter-filled umklapp 
scattering becomes always dominant at high frequency 



when the parameter is close to the boundary. There- 
fore the power-law dependence of the optical conductiv- 
ity <j(lo) oc uj 16K p- 5 could be compatible with the charge 
gap, which originates in the Mott insulator and/or the 
CO insulator. 

In conclusion, we found two kinds of insulator, the 
Mott insulator and the CO insulator, for the ground state 
of a one-dimensional quarter- filled system with dimeriza- 
tion, where the half-filled umklapp scattering induced 
by the dimerization competes with the quarter-filled 
umklapp scattering for large nearest-neighbor Coulomb 
interaction V(> U/A). In the Mott insulating phase 
close to the boundary, the half-filled umklapp scatter- 
ing is dominant at low T but the quarter-filled umklapp 
scattering becomes dominant at high T, and then the 
resistivity could exhibit a crossover from a quarter-filled 
behavior to a half-filled one as decreasing temperature. 
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